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LIST  OF  SYMBOLS 

Piston  area 
Primary  pressures 

Secondary  (difference  frequency)  pressure 
1  +  B/2A  nonlinearity  number  a  3.6  for  water 
Normal  density 
Normal  sound  speed 
Secondary  wavenumber 
Geometric  mean  primary  wavenumber 
Geometric  mean  primary  Rayleigh  distance 
Range  (vector) 

Range  (scalar) 

Angle  to  field  point 

“l  +  a2  =  Absorption  (nepers/m) 

Integration  range  variables 
Scaled  range 
Absorption  number 
Nearfield  absorption  number 

Farfield  absorption  number 
Downshift  ratio 
Westervelt  source  level 

PL/W  where  PL  is  the  apparent  source  level  at  L 

u  +  iv  complex  variable 

/4«/k'  =  Westervelt  angle 

0/0 
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ABSTRACT 


It  is  often  necessary  to  measure  and  operate  parametric  radiators 
under  conditions  where  finite  range  effects  are  important.  A  nearfield 
model  is  presented  in  which  the  source  level  is  obtained  by  assuming 
the  primary  waves  to  be  plane-cylindrical  or  diverging-conical  depending 
upon  whether  the  observation  range  is  less  or  greater  than  the  Rayleigh 
distance.  Asymptotic  beam  patterns  are  obtained  from  the  solution  of 
an  exponentially  tapered  line  array.  Nonlinear  attenuation  is  not  consider¬ 
ed. 
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This  memorandum  was  prepared  under  Project  A-614-15,  Subproject  11-121- 
706  Task  17442  "Parametric  Echo  Ranging  Systems",  Principal  Investigator 
W.  L.  Konrad,  Sponsoring  activity  NSSC  SEA  06H1  J.  Neely. 
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INTRODUCTION 

In  the  Westervelt^  model  of  the  parametric  radiator,  an  asymptotic 
approximation  was  made  to  obtain  the  source  level  and  beam  pattern  at  very 
long  ranges.  In  many  practical  cases  the  conditions  for  this  model  are  not 
met,  therefore  a  finite-range  model  is  essential. 

2-  7 

Several  methods  '  have  been  employed  to  obtain  the  field  at 
finite  range  using  various  simplifying  assumptions.  In  the  present  method 
the  asymptotic  beam  pattern  for  points  sufficiently  far  outside  the  primary 
beam  is  obtained  from  the  solution  for  an  infinite  exponentially-tapered 
line  array.  The  source  levels  are  obtained  by  assuming  the  primary  waves 
to  be  plane-cylindrical  or  diverging-conical  depending  on  whether  the 
observation  point  is  less  or  greater  than  the  Rayleigh  distance.  The  other 
assumptions  are  that  the  absorption  is  linear  and  small  compared  to  the 
secondary  wave  number. 


Beam  Patterns  of  a  Line  Array 

1  For  a  parametric  array  the  source  strength  per  unit  length  is  given 


by 


q(x,t)  =  0$  d/dt  (P1P2)/ooV 


(1) 


where  Pn  =  PQn  exp  [-x(«n  +  ikn)  +  iu>nt],  n  =  1,2 

The  difference  frequency  pressure  for  small  S  becomes 

00 

P(r,t)  =  -pq  j  dx  q*  (x,t  -  Jr”  -  x|/c0)rtnr  |‘r  -  x  |  (2) 

0 

where  q'  indicates  the  time  derivative 

For  rU«  Eq  (2)  gives  the  Westervelt  result 

rP(0)  =  PP1P2k2S/8’r  pqc02|“+  ik  sin2(0/2fj  (3) 

and  W  =  rP(O)  is  the  Westervelt  source  level. 

(The  term  exp  (i^t  -ikL)  is  omitted  throughout.) 

For  finite  range  Eq.  (2)  can  be  written  in  terms  of  the  dimensionless 
variable  Y  =  PL/W  according  to  the  geometry  of  Fig.  1  as 
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00 


/ 

-L 


dx 


exp  [-2*x  -  1k(v^~ h  a2  +  x)] 


Vx2  +  a2 


(4) 


Eq.  (4)  can  be  rewritten  as  a  contour  Integral 


.  where 


Y  =  2«L  exp(2<*L)  f  d?-  -exP(.~z) 

J  JFTp 

o 

z  =  2«x  +  ik(v^  +  a2'  +  x) 
20«-2«L  +  ika2/2L 
B  =  2a  (“2  +  ik«01/2 


(5) 


and  a  typical  contour  is  shown  as  the  upper  curve  of  figure  2. 

Since  k/  a  »i  is  the  domain  of  interest  the«  term  in  B  will  be 
neglected.  Further,  for  small  angles  it  Is  convenient  to  normalize 
0  «  a/L  to  the  Westervelt  angle  ft.,  hence, 

W 


z  =  +  iv  M  a*L(-l  +  ie2) 

ooo 

B2«  4i(2«L0)2  =  4iu0v0 

where  0  =  0/0  . 

w 

Since  the  path  of  integration  becomes  rapidly  oscillating  for  large 
z,  the  lower  path  of  Fig.  2  where  v  =  v  may  be  substituted  to  give  the 
same  result. 

Changing  variables  to  z  =  t  -  uQ  +  v  .  Eq.  (5)  becomes 

00 

Y  =  2«L  exp(-ivQ)  f  dt.e*P(-t?  .  (6) 

0  'A*  -  u0)2  -  v„2  *  Z,vo(t  +  uo>’ 
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Eq.  (6)  was  evaluated  numerically  and  the  results  are  shown  in  Figs.  3  and 
4.  Note  that  for  the  line  source  (S->0)  there  is  a  logarithmic  singularity 
of  strength  exp(-2<*L)  as  6-*0.  Only  for  values  of  2o-L>10  do  patterns  fall 
within  1  dB  of  the  Westervelt  asymptote.  For  smaller  values  of  2«L  the 
patterns  tend  to  become  narrower,  approaching  a  minimum  at  8  =  1  and 
Increasing  thereafter0. 

The  line  array  patterns  are  asymptotic  patterns  of  a  parametric 
radiator  at  points  sufficiently  far  outside  the  primary  beams.  For  experi¬ 
mental  comparison  the  axial  level  must  be  known.  In  order  to  determine 
the  axial  level  it  is  necessary  to  account  for  the  finite  cross-section 
specifically. 


CYLINDRICAL  SOURCE 

For  a  plane  wave  cylindrical  source  of  radius  a  and  the  field  point 
on  the  axis  we  have 


.  PPlP2k2exp(2«L) 


4  "■  ?  c  ‘ 
0  0 


®  a 

/dx  / 


-L  0 

Integrating  over  a  we  have 
/JP^k 


da  2ira  exp[-2<*x  -  ik(v>c2  +  a2'  +  x]  ,7% 

v^77 


r  o 


p  = 


1  2PnCn 
0  0 


2  exp(2«L) 


/ 


00 

dx  exp(-2«*x)  +  J  dx  exp(-2«x  -Jikx) 
-L  0 


00 

j '  dx  exp(-2«x  +  ik(Vx2  +  a2'  +  x) 


-L 


(8) 


The  second  integral  is  omitted  since  k  »«  and  its  contribution 
negligible.  Taking  z~=  2<*x  +  ik(\/xz  +  a2'  +  x)  in  the  third  integral  and 
neglecting  terms  in  we  obtain 

4»x  =  z  -  vz2  +  B2' 


from  which 


4<*dx  = 


dz  £l 
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Substitution  of  the  modified  contour  z  =  t  -  u  +  1v  in  the  third 
Integral  gives  finally  w 


Y  -  1 

Y°“  TT 
o 


-  \  exp(-ivQ)  J  dt  exp(-t) 
0 


t  -  u„  +  ivrt 
_ o _ o 

'/(t  -  “o>2  -  »/ 


+  2iv0(t  + 


(9) 


where  vQ  =  ka^/2L  2  oc  Rq  k/ko  Uq  .  .  ’ 

Eq.  (9)  was  evaluated  as  a  discontinuous  function  at  t  =  u  ,  since 
the  integral  converges  more  rapidly  than  by  integrating  the  terms  separately. 


The  results  shown  in  Fig.  5  give  the  apparent  source  level  as  a 
function  of  scaled  distance  referred  to  the  Westervelt  source  level  at 
infinity.  The  parameter  is  the  nearfield  absorption  number.  Oscillations 
occur'  as  the  source  level  vs  range  changes  from  square-law  to  a  slower 
rate.  For  all  values  of  the  parameter  the  asymptotic  approach  for  the 
Westervelt  limit  is  from  above.  This  effect  has  a  simple  geometric  explana¬ 
tion,  i.e.,  the  distance  is  measured  from  the  origin  rather  than  the  "center 
of  gravity"  of  the  array,  hence  apparent  source  levels  tend  to  be  higher. 

Note  that  the  effects  of  absorption  are  negligible  for  2«L  <  0.1. 

From  2«L  =0.1  to  1  absorption  begins  to  have  an  effect  and  for  2*L  >  1 
The  curves  approach  the  Westervelt  limit.  The  resulting  curves  are  asymp¬ 
totic  for  L«R  . 


CONICAL  SOURCE 


For  l»R  consider  a  conical  source  of  spherically  diverging 
primary  waves  and  field  point  on  the  axis.  Eq.  (7)  becomes 


PlP2kV 


2p  c 
0  0 


0  ® 

J *  d9sin6  J"  dr 

0 


■[2<*r  +  ik 


(Vr2  +  L‘ 


2Lrcos9  +  r  -  L)] 


Vr2  +  L2  -  2Lrcos0 


(10) 
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Integrating  over  6  we  obtain 


P  = 


WoV 


1  2  *nCn^L 
0  0 


r 


dr  exp[-2arr  -  2ik(r  -  L)] 
r 


% 


A  7 

exp[-2ar  -  ik(^r  +  L  -  2Lrcos0_  +  r  -  L ) ] 
dr  - - - ^ -  (11) 

0 

The  second  integral  is  again  omitted.  Substituting  z  for  x  in  the  third 
integral  as  in  the  cylindrical  case  and  taking  the  new  contour  we  have 
after  some  manipulation 

uo 


f  exp(-t) 

-J  « -  -  . 

0  1  +  uo  +  iv0  -J  (t  -  u0)2  -  vo2  +  2ivo(t  +  \]  J 

where  vD  -  k  L  0^/2  at  1^/(2  ^Rq  ko/k)  and  Res.  is  the  residue  at  zD  . 

Although  each  integral  has  a  pole  at  t  =  0,  the  sum  is  well-behaved. 
Near  zero  the  last  integrand  of  Eq.  (11)  behaves  as  1/z,  where  z'  = 

(1  +  iv  /u  )t  and  it  is  necessary  to  evaluate  the  residue  between  the  two 
paths.  °Taxe  the  contour  from  z '  =  (1  +  i v  /u  )«  to  z1  =  (1  +  i0)$ 
and  let  «  -*  0,  which  gives  0  0 

/  fr'  *  -  7  In  [1  +  (v0/u0)2]  -  i  tan-^v^) 

C 


r° 

Res.  ♦  /  dt 
0 


/  (t  -  uj2  -  V  2  +  2i v  (t  +  uj 


Using  the  residue  as  the  initial  value,  Eq.  (12)  was  integrated  as 
before  and  the  results  are  shown  in  Fig.  6.  The  parameter  in  this  case 
is  the  farfield  absorption  number.  For  values  greater  than  1,  the  approach 
to  the  asymptotic  is  again  from  above.  For  values  less  than  1,  the  asymp¬ 
totic  source  levels  fall  below  the  Westervelt  asymptotic  and  approach  the 
asymptotic  from  below.  This  is  the  farfield  diffraction  effect  where  the 
generation  becomes  spherical  and  the  source  level  approaches  the  exponential 
integral  rather  than  the  Westervelt  asymptote8. 
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Figure  7  showSgthe  predicted  source  levels  vs  range  compared  to 
experimental  results5.  The  downshift  ratio  was  7  and  the  value  of 
2«R  was  0.0175.  The  cylindrical  and  conical  curves  intersect  near  the 
RayTeigh  distance.  Above  this  point  agreement  with  the  conical  curve 
is  excellent  while  below  the  deviations  are  significant. 

Figure  8  shows  a  comparison  of  experimental  and  predicted  source 
levels  vs  range  for  a  downshift  of  50  with  2*R  =  0.17.  Experimental 
values  are  approximately  3  dB  low  near  the  RayTeigh  distance. 

Figure  9  presents  beam  patterns  of  the  second  example  made  at  two 
different  distances:  2«L  =  0.25  and  1.  In  the  first  case  the  experi¬ 
mental,  rather  than  predicted  axial  value  was  used. 

CONCLUSIONS 

Predictions  based  on  the  nearfield  model  appear  to  be  in  fair 
agreement  with  experiment,  however,  further  comparisons  must  be  made  to 
determine  all  the  conditions  for  validity  since  errors  can  be  large 
near  the  Rayleigh  distance.  The  particular  problem  of  nonlinear  atten¬ 
uation  must  also  be  included  if  the  model  is  useful  in  a  practical  sense. 
In  its  present  form,  however,  it  serves  to  illustrate  the  various  finite 
range  phenomenon  of  an  absorption  limited  parametric  radiator. 
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FIGURE  CAPTIONS 
Geometry  of  Line  Source 

Contour  Integration  Paths.  Branch  points  lie  outside  both  paths 
and  are  not  involved. 


Beam  Patterns.  9  is  the  angle  referred  to  the  Westervelt 
angle  and  the  parameter  is  the  scaled  distance  2*L  =  10  to  1 

Beam  patterns  2“L  =  1  to  .01 

Source  level  relative  to  Westervelt  vs  scaled  distance. 

The  parameter  is  the  nearfield  absorption  number  (plane 
wave  model ) . 

Source  level  relative  to  Westervelt  vs  scaled  distance. 

The  parameter  is  the  farfield  absorption  number  (spherical 
wave  model ). 

Comparison  of  predicted  and  experimental  values  of 
apparent  source  level  vs  range  2«Rq  =  0.0175,  kQ/k  =  7 

Comparison  of  predicted  and  experimental  values  of 
apparent  source  level  vs  range  2®r0=  0.17, ko/x  =  50. 

Comparison  of  predicted  and  experimental  beam  patterns 
2«L  =  0.25  and  1,  kQ/k  =  50. 
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Fig .  2  -  Contour  Integration  Paths 


Fig.  1  -  Geometry  of  Line  Source 
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Source  level  relative  to  Westervelt  ■ 
distance.  The  parameter  is  the  near 
absorption  number  (plane  wave  model) 


Corrparisor.  cf  predicted  and  experi-uental  values 
of  apparent  source  level  vs.  range  2  = 


